function out = j_r( r )
%J Summary of this function goes here
%   Detailed explanation goes here

global q;
global u;
global inf_jt;
global h_jt;

tmp = -q ./ u;

out = tmp .* trpzInt(@(t)int_jt(t, r), -inf_jt, inf_jt, h_jt);


end

